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ABSTRACT 

We address theoretically adiabatic regime of charge transport for a model of two tunnel-coupled quantum dots 
connected in series. The energy levels of the two dots are harmonically modulated by an external potential with a 
constant phase shift between the two. Motivated by recent experiments with surface-acoustic-wave excitation, we 
consider two situations: (a) pure pumping in the absence of external voltage (also at finite temperature), and (b) 
adiabatic modulation of the current driven by large external bias. In both cases we derive results consistent with 
published experimental data. For the case (b) we explicitly derive the adiabatic limit of Tien-Gordon formula 
for photon-assisted tunneling and compare it to the outcome of simple conductance modulation. A tutorial for 
adiabatic pumping current calculations with the Green function formalism is included. 
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1. INTRODUCTION 

Alternating external fields applied locally to a confined electronic structure can result in a measurable time- 
averaged current across the device. Generating dc current by microwave excitation at frequencies of several GHz 
and less is often referred to as pumping. Particular appeal of pumping both as a research tool and a potential 
source of applications for nanoelectronic systems is that low-frequency fields can be confined to wave-guides and 
wires and thus conveniently delivered to samples at a few Kelvin or lower temperatures. This opens the way to 
explore pumping in situations when quantum coherence and single-electron charging effects matter. 

Particular motivation for this work has been provided by two recent experiments^'^ in which electrons in 
quantum dots have been subjected to an alternating piezoelectric potential of a running surface acoustic wave 
(SAW). The experiment of Buitelaar et al} focused on pure pumping current with no external dc voltage 
bias. Experimental results have been reported to be in a good agreement with the theory of adiabatic quantum 
pumping'^ appHed to a simple two-level model (see definition and discussion in Secs.[2]and[3]of the present paper). 
The experiment of Naber et al.^ has been designed to measure the influence of an oscillating SAW potential on 
directed current driven by a constant external voltage but not necessarily by the SAW-induced potential itself. 
The results were found to be in a good agreement with Tien-Gordon formula"*' ^ of photon-assisted tunneling. 

In this paper we show that the results of these two experiments are consistent within one and the same model: 
in the large bias limit, adiabatic conductance modulation for the two-level system of Ref.^ reproduces the major 
features of the current traces reported in Ref.^ (except for the fine structure due to non-adiabaticity). We 
derive analytically the adiabatic limit of Tien-Gordon formula and reproduce the criterion® for the validity of an 
adiabatic approximation: inverse life-time of the discrete charge states, P/ft, must be larger than the modulation 
frequency lo. 

This paper has two major components. After introducing the model in Sec. [21 we discuss in detail the 
adiabatic pumping current and its various limits (Sec. [3]). These results have been used directly for interpretation 
of experimental data on pumping in carbon nanotube quantum dots.* The discussion in Sees. [2] and [3] is more 
of a tutorial style. In the second part of the paper. Sec. [H we turn to the limit of large bias in which adiabatic 
conductance modulation dominates over pure pumping. We compare predictions of the adiabatic time-domain 
expansion^ (essentially, time-average of the Landauer formula) to the low-frequency limit of the photon-assisted 
tunneling through a double dot.*'^ This comparison uncovers precise agreement between the two theories (up 
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to an overall integer factor) and establishes a rigorous adiabaticity criterion. Comparison to the experimental 
data of Naber et al identifies the adiabatic features in the measured current line-shape. A brief summary and 
an outlook at open challenges in Sec. [5] concludes the paper. 



2. THE MODEL AND ITS STATIC PROPERTIES 



The Hamiltonian Ti. describes two energy levels, ei and £2, with off-diagonal tunneling coupling A/2, connected 
to two external reservoirs: level ei to the left (L), level £2 to the right (i?). In the second quantized form, 



n=J2 + (4^2 + 4^1) (A/2) 
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The boundary conditions imposed onto Ti. are those of equilibrium reservoirs at electrochemical potentials fiL and 
(for the left and the right lead correspondingly) and at equal temperatures T. The difference fi^ ^ fJ'R = eVbias 
specifies the external dc bias, while the average (/il + I^R.)/'^ = n ~ sets the reference level for energy. 

We shall take Vka in Eq- (HI) to be independent of energy (the limit of a structureless, wide band ). In this case 
the only parameter characterizing the leads is the golden rule half-width = 27r J^k l^fea P'^(Mq ~^fea)- If no lead 
index is specified, the couplings will be assumed symmetric, — Fj? = F. When discussing near-equillibrium 
properties (linear conductance, abiabatic pumping) we shall assume vanishing bias conditions, Vbias 0+. 

The reason we omit physical spin index in Eq. ^ is that the dots are assumed to be in the Coulomb blockade 
regime, so that double occupancy of a single dot in energetically prohibited. This assumptions of "spinless & 
non-interacting" electrons can be viewed as a fermionic representation of the N N + 1 charge state transition 
in an otherwise inert (that is, devoid of any internal dynamics) Coulomb-blockaded quantum dot. For a single 
dot this representation results in the well-known single resonant level (Breit-Wigner) approximation to Coulomb 
blockade.® Applicability of a non-interacting Hamiltonian like ([T]) to a pair of tunnel-coupled quantum dots is 
thoroughly discussed in Sec. V of a review^ by van der Wiel et al. While incomplete in its treatment of mutual 
capacitance and correlation effects, the non-interacting charge-carrier approach adopted here allows for a very 
detailed investigation of time-dependent effects because of the its fundamentally single-particle nature. 

Some aspects of non-adiabatic pumping for the system defined by Eq. Ilj have been considered recently'^^ 
using a Floquet formulation for time-periodic non-equilibrium Green functions. Our model is also a special (and 
previously unexplored) case of the a tight-binding model for SAW-induced adiabatic pumping as introduced and 
studied in Ref.^-^Q 

It is convenient^^ to define the retarded Green function matrix, G, for the double-dot region in the following 
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way,^^ 
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The transmission probability from left to right. 
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determines the linear conductance via (the finite-temperature) Landauer formula, Q — {e^/h) J dE{~d f / dE)T {E) . 
Here and below fa{E) = [l + e^-^^^")/'^^-'"] ^ is the Fermi distribution function. 
For our system the transmission probability is [see Eqs. ^ and ^] 

A2 TlTr 



r{E) 



AD{E) ' 



(5) 



*In the notation of Ref.,^'^ one has to consider N — 2, take the limit of J — > 00 with Ja/ J ~ const and identify Jd 
with A/2. 
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Figure 1. (Color online) Schematic plot of the stability diagram^ for our model as function of £i (in units of A). Dashed 
(yellow) lines show the eigenenergies Ei for no coupling to the leads (F = 0). The tunneling splitting A is the maximum 
of El — E2. The grayscale plot is the conductance Q for F/A — 0.2, ksT j ^ — 0.1, darker shade corresponds to higher 
values of Q. The number in the parenthesis characterize the charge state for each of the dot. For F/A exceeding 1 the 
two dark blobs (triple points) merge into a quadruple point (not shown). 



were the denominator DiE) = | detG'(i?)| ^ is given explicitly by 

D{E) = {[E ~ e{\[E - £2] - AV4)2 + {E ~ eifvl/A + {E ~ £2)'ri/4 + TlVrI\''/% + riryi6 (6) 
^{[E-E^f +T^/A){[E-E2f +V'I'^) ioYTL = TR. (7) 



The eigenenergies Ei and E2, £'1,2 = (ei+£2±-\/(^i ~ £2)^ + A2)/2, can be understood as "molecular orbital" 
states of the double dot system. For F < A the resonance lines of i?i,2 = are resolved in a two-dimensional plot 
of T{E = /Lt) as a function of £1 and £2, see Fig.[TJ Such diagrams are often referred to as stability diagrams^ for 
systems of coupled quantum dots since they identify the regions of different equilibrium occupation numbers for 
the dots. Anti-crossing behavior seen in Fig. [T] is a hallmark of a tunnel-coupled double dot system and serves 
as a conceptual basis for charge-qubit approaches to quantum information processing. 

3. ADIABATIC PUMPING 

We envision pumping by modulating £1 and £2 with the help of a piezoelectric potential of a propagating SAW 
(or, perhaps, by direct gating). This requires calculating the response of the system to time-variation of the 
double-dot parameterfQ. One of the simplest scenarios^ is a phase-shifted harmonic modulation: 

ei{t)=eo + 6/2 + Pcosujt, (8) 
£2 {t) ^eo-S/2 + P coa{ujt + ip) . (9) 

where 5 is a constant detuning between the levels and P is the amplitude of the external potential. However, 
in this section we shall focus on general properties of the adiabatic pumping current which are independent of a 
particular choice of Si,2{t). 

Let us apply the general theory of adiabatic transport in coherent strctures'' to Hamiltonian ([T]) using a Green 
function formalism. Instantaneous adiabatic current XP"'"P(i) entering lead a from the nanostructure is^^ 



_p.,^ _± f I ~d[hiE) + ME)] 

" ^ ' " 27r / 2 dE 



G(i?)f„ G^E) jH" 



dE . (10) 



^For simplicity, the effect of the modulating field on the tunnel couplings is assumed to be negligible. 
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Figure 2. (Color online) The response function R{ei,e2) allows for determination of the charge pumped per period for an 
arbitrary contour. Blue dashed line shows an example of a pumping contour that produces the average current close to 
one electron per cycle. The model parameters are F/A = 0.2, fcsT/A = 0.1; ei,2 are in units of —A. For these values, 
the maximum possible value for the pumped charge per period is Qmax = 0.981e. 



The charge pumped by the system during one cycle of the periodic modulation with a circular frequency u> is 

-df{E) 



Jo 



'T-pump 



{t)] dt = e 



dE 



[Pi {ei,e2,E)dei + P2{ei,e2,E) de2] dE , 



where 



Pi{si,e2,E) 
P2{ei,e2,E) 



IQttD{E) 



(11) 

(12) 
(13) 



1QttD{E) 

Since we have used a symmetrized current expression in Eq. pT|) . Pi ^ — F2 under 1 2, L ^ i?. 

As usual for two-parameter pumps, ^ it is convenient for the subsequent analysis to transform the contour 
integral in Eq. pT|) into an surface integral using Green's theorem. We write Q as a double integral over the 
area enclosed by the contour: 



Q ^ ±e JJ R{ei,e2) deide2 , R= J 



f -dfiE) 




dPi 


J dE 


dei 


de2 



dE. 



(14) 



Volume under the two-dimensional plot of R{si, £2) above the pumping contour's image in the R — plane gives 
the value Q. The sign of Q is determined by the direction of the contour [plus in Eq. for counterclockwise 
direction] . This way of visualizing the pumping response is illustrated in Fig. [21 

Restricting the discussion to F^ = Tji allows for some further analytical progress. Performing the integration 
over the energy E, we obtain explicitly 

A2 Im [S'iEi] - 5'(^2)] A^F [{E2 - Ei) Rc{S"iEi) - S"iE2)} + rim{S"iEi) + S"iE2)}] 



2tt{Ei - E2) 



An{Ei-E2Y[{Ei-E2Y+V'] 



(15) 



where S{E) ee 'I' [1/2 + (F + 2iE) / {-i-KkBT)] and ^' is the digamma function. 

Let us consider also the maximal possible charge per period (Qmax) which is obtained when a large pumping 
contour completely encompasses the positive part of R. For an arbitrary ratio T /kgT we find: 



A2 {2(a;2 + F2) ImS'(a;/2) -f xF [a:ReS"(x/2) - FImS"(j;/2)]} 



/A 7ra;2(a;2-HF2)Va;2- A2 

Specific limits of R and Qmax are illustrated below. 
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Figure 3. (Color online) The maximal pumped charge per period in the limit of T /ksT (left) and ksT /T (right) 
as a function of fcsT and F respectively. For three representative ratios F/A = 0.3, 1 and 3 (left) and ksT / A — 0.1, 0.4 
and 2 (right) the behavior of the response function R is shown (absolute scale of R varies from inset to inset). The values of 
Qmax for the insets are (0.976, 0.646, and 0.146) and (0.994, 0.671, 0.183) in the tunneling- and the temperature-dominated 
regimes respectively. 



3.1 Tunneling-dominated limit: fc^T <C T 
The function R{ei,ei) for ksT <C T snTiplifies to 

i?(T ^ 0) = - 

The maximal possible charge (|T6l) becomes 

Qmax(r 0)/e = 1 - 



r^A^(£l+£2) 

87rL>2(0) 



Vr2 + A2 



1 - (r/A)3 r < A 
3AV(2r2) r > A 



(17) 



(18) 



This is illustrated in the left panel of Fig. [31 We see that charge quantization (meaning Q — > e) is possible when 
r <C A. In this limit the two resonance lines in Fig. [T] are well-defined, and pumped charge quantization proceeds 
via the loading- imloading mechanism. 

3.2 Temperature-dominated limit: ksT ^ T 

In the opposite extreme of temperature-broadened resonance lines we have 



(prime denotes energy derivative). The corresponding maximal charge is 

r tanh[.A/(4fc.r)] ^^ ri-C^(e-A/('--^)), fc.T « A 
^ ' " ^ ~ |^A/(8fcsr), fcBT»A 



(19) 



(20) 



This limit of the model approximates the conditions of the seminal experimental work of Pothier et al}^ Note 
that the "exponential" accuracy is possible only as long as F ^ fc^T so that co-tunneling and other higher-order 
processes remain below the level of thermal fluctuations. 
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4. CONDUCTANCE MODULATION AT LARGE BIAS 

4.1 Adiabatic pumping in the presence of external dc bias 

Let us turn to a situation when a constant bias voltage Vbias is applied together with a slow modulation of 
the system parameters. General formalism for this scenario has been considered in Ref.^ The instantaneous 
adiabatic current X{t) from left to right has been found to contain two additional components on top of the pure 
pumping contribution (jlOp : 

I{t) ^ XP"'"P(t) . (21) 

In our notation, the results of Entin-Wohlnian et al7 for the time-average current read 

^^^^C^^^^l J [ME)-fRiE)]nE)dE, (22) 

^ ' 2tt Jo h J 2 dE ^ ' dt ^ ^ 

Here (pr is the overall transmission (Fricdel) phase. 

Interpretation of Eq. ([221) is straight-forward'' — it is the time-average of the current given by the finite- 
bias Landauer formula. Adiabatic perturbation of the system manifests itself in Eq. (|22p only as a parametric 
modulation of the transmission function T(E). The last remaining term, Eq. (|23p . is harder to interpret. It is 
non-zero only if both the bias Vbias 7^ and the adiabatic modulation are present at the same time (thus the 
notation X""'^). 

We observe that both IP"'»p [Eq. ([TO])] and X'"''' [Eq. ^] contain the Fermi functions fa{E) only in the 
form of the respective energy derivatives. If eVbias 3> fc^T and the bias window E G [iJiu, hl] is sufficiently wide 
to cover the whole energy range where the transmission function T(E) is appreciable, then the conductance- 
modulation term, I^^'^, dominates in the total current (HI]). In the remaining part of this section we are going 
to explore the large bias limit eVbias ^ \^i,2{t) ~ fJ'\,kBT, Tl^h for the model defined in Eq. ([!]) and subjected to 
a harmonic modulation ([S]) and 

4.2 Large-bias current through a modulated double-dot system 

The average dc current in the large bias limit, /^^lab _ jbias^^'j^ jg 

.adiab ^ £ frmdF ^ ^ r^'^ AHTL+Tn)dt 

hj ^ ' h Jo (ri + rfl)2(i + AV[rLrK]) + 4[ei(t)-e2(i)]2 ■ ^ ^ 

To derive the above equation we have performed the energy integration in Eq. (|22p using Eq. ^ and extending 
the bias window to ±00. The result for the time-averaged directed current ([24]) is of the form 



J "'Tw^h^- '''' 

with Peff = 2P sin (^/2 and some positive constant A. Thus the comparison to an experiment can be made without 
the detailed knowledge of the tunnel couplings. Only the three scale factors — for the current /, the detuning S 
and the modulation amplitude Peft — have to be adjusted. We compare Eq. (|25|) with the experimental results 
of Naber et al in Fig. U) Our result quantifies the qualitative adiabatic arguments of Naber et al. concerning the 
increased time at resonance near the turning points of harmonic modulation.^ 
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Figure 4. (Color online) Comparing the large bias limit of the adiabatic theory [Eq. (|25p . blue continuous lines] to the 
experimental data^ (black circles) on a double quantum dot system formed electrostatically along an etched channel in 
a GaAs/AlGaAs heterostructure. The units are arbitrary (for quantitative details of the experiment, see the original 
work by Naber et al. ). The horizontal axis is the difference between the gate voltages on the two dots (proportional to 
S), the vertical axis is the current measured in the presence of a fixed bias voltage (same for all graphs) and a certain 
relative amplitude P^s due to applied SAW power. The latter is zero for the lowest trace and increases from bottom to 
top (individual traces are shifted along the /-axis for clarify). All graphs share the same scales for 5 and for Pcff, but 
the overall magnitude of the current is fitted separately for each P. Note that the ripples between the broad peak are 
non-adiabatic features^ that are adequately described by Tien-Gordon formula,^ see Eq. (^5)1 . 



4.3 Comparison to the non-adiabatic theory of Stoof and Nazarov 

Large bias limit for transport through two capacitatively coupled quantum dots connected in series has been con- 
sidered theoretically by Stoof and Nazarov.^ In the limit of a very weak inter-dot tunneling, A <C F^, r^, hu), PgS, 
zero temperature and strong asymmetry ,Tl<S^Tr, they find the following expression for the time-average cur- 
renlll 

Ah ^ r2„/4 + (^ - nhwY ■ ^ ' 

Here J„ is the n-th order Bessel function of the first kind. A similar formula is known from the theory of photon- 
assisted tunneling between two superconductors due to Tien and Gorgon,'* therefore Eq. (j26p is sometimes 
referred to by their name.^ It has been successfully applied to several experiments on microwave excitation of 
double quantum dots,^ as well as to the measurements of Naber et al? that we used in Fig. [4] 

We would like to compare Eq. (|25p to — s- limit of Eq. (j26p . This limit is singular and requires some care 
in implementation. To this end, we note that Eq. (j26p represents a convolution of the Bessel function with a 
Lorentzian, see Fig. [5l Typical width of this Lorentzian along a (dimensionless) n-axis is V {fno) while the 
period of oscillations of is of order one. Evidently, in the limit Vn ^ Tilo the fine oscillatory structure of the 
current as a function of 5 will be washed out. This is the validity condition for the adiabatic limit to Eq. (|26p 
that we are in a position to derive. 

•'"This is equation (21) of Ref.^ in our notation. 
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Figure 5. (Color online) Derivation of the adiabatic limit to the theory of Stoof and Nazarov [Tien-Gordon formula, 
Eq. ([Mil]. 



The Bessel function Jn{z) oscillates for z < n asd 



Jn{z) 



COS [z9{n/ z)] + sin [z9{n/ z)] 



^2)1/4 



for \z\ < n , where 6{S,) = — £, arccos^ . 



(27) 



Equation (P7)) contains the analytic form for the envelope of oscillations shown in Fig. [5l Replacing the sum in 
Eq. (|26p with an integral, using Eq. (j27p and dropping the rapidly oscillating term gives 



I'^iu; 0) 



e 



-1/2 



e 

Ah 



A^TrcIt 



(28) 



/_! 72/4 +{S- CPeff )2 4h Jo ry4 +{S- Peff sin r)2 ■ 

By derivation, Eq. (|28p approximates well the full non-adiabatic result (pS]) if < huj. 

If we now expand the adiabatic result to the leading order in A and take the limit F^ <C Fjj then 
Eq. is recovered up to a factor of 4, 



(29) 



To a reasonable degree of confidence we have excluded the possibility that the factor of four arises due a 
notational or/and algebraic mistake. Of course, we should bare in mind that the Hamiltonian of Stoof and 
Nazarov is different form ours: they consider the case of both strong intra-dot and inter-dot Coulomb repulsion. 
Only three states of an isolated double dot system are possible in their case while in ours there is the forth state 
— the doubly occupied configuration (1, 1). 

5. CONCLUSIONS 

Steady progress in experimental manipulation of nanoscale systems such as quantum dots and wires has recently 
reached a level enabling direct comparison between theory and experiment for pumping in these devices. Many 
observed features of the pumping currents can be explained within a framework of single-electron models with 
time-dependent energy levels. In this paper we have looked at two examples of application for this framework 
to recent experiments. It is quite satisfactory to find agreement between different limits of different devices 
that comes within a single model. It is conceivable that in the future more refined experiments will allow 
for quantitative characterization of the pumping current magnitude beyond the accuracy of a simple single- 
electron picture. Clear separation between the effects of discrete charging, quantum interference and many-body 
correlations in pumping through nanostructure seems posees future challenges both for theory and experiment. 



^The highly accurate approximation (|27|l is derived by taking n = and applying the stationary phase expansion in 
to an integral representation of the Bessel function, Jn{z) — {2-k)~^ J^^ exp{iz sin (p — irvp) dip. 
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